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1. Introduction 


Contract NAS8‘-34768 callad for the developaent of elgorithas for 
correcting certain aeaaureaent art if acta present in the aeasureaents aade 
during the 1981 flight prograa. This problea was addressed in Technical 
Report no. 1, issued in April I982y and subsequent progress reports. The 
purpose of this supplement is to suaaarize the recoomended algorithm for the 
correction of higher-order data errors. 

II. The source of higher-order errors 

As noted in Technical Report no. 1, there are several sources of error 
in addition to the drift-angle error treated in that report. These include 
the non-uniform attitude-measureswnt delay (due to the non-uniform scanning 
intervals), possible laser drift, and laser scanning above or below the 
horizontal. Due to their random natures these errors cannot be corrected by 
deterministic means. 

III. Suggested correction algorithm 

Fortunately these higher-order errors are essentially independent from 
scan to scan. As such the error spectrum they produce is white, and error 
reduction can be obtained by recourse to spectral filtering. 

Where the errors are constant with range, as in the case of most of the 
instrumental errors, a continuity argument can be used to justify a simple 
correction algorithm. Consider a time series of mean radial-velocity 
measurements. Each item in the time series is an average ove. aboitt 10 km of 
range. The mean-square difference between adjacent items in the series will 
arise from that portion of the two-dimensional velocity spectrum which is not 
averaged out by the 10-km one-dimensional integration. That is, the major 
part of the variance of the true time series will be due to spatial 
variations in the wind field with scale sizes on the order of several km. Due 
to the integration the spectrum of this true time series will be concentrated 
at low frequencies. 

The observed time series will contain in addition variance due to the 
errors mentioned above. Since the errors are essentially independent from 
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«c«n to scan, tha additional varianca in tha apactrun vill ba whita. Sinca 
theaa two contrl'witiona to tha apactrua differ greatly in apactral fora, 
filtering tachniquaa can aubatantially reduce tha white coaponant. A 
transveraal low^paaa filter la uaed to eatiaate tha true coaponent of tha 
time aeriea* Tha paaaband ia choaen to be aufficiently broad to paaa all 
aignificant frequency coaponanta preaent in tha integrated velocity 
aeaaureaMnta. Thia ia eaaily accoapliahed aince the frequency range of tha 
true tiaa aeriea haa been greatly reatricted by integration along the lidar 
beaa. 

Many typea of filtera could be uaed to implement thia algorithm. Tha 
following auggeation makea uae of a leaat-aquarea fit of a cubic function to 
the data, to allow for alow changes in the mean radial velocity. Forward and 
aft scan data sets are considered separately. Radial velocities are averaged 
over the longest range interval over which the measurements are highly 
reliable (by signal'-to-noiae ratio criteria). Each point in the resulting 
time series is treated independently: a correction for what point is derived 
from adjacent uncorrected averages, and that correction is applied to 
measurements at each range g4te. 

The corrected swan radial velocity for a g:ven point can be derived as a 
simple weighted mean of the mean apparent velocities in the neighborhood of 
that point. Using an odd number of points (2M-*-l) centered on the point in 
question, let the mean radial velocities be Vi, i"-M to M. For cubic or 
quadratic approximation (they give the same result for symnetrical 
filtering), the corrected mean velocity V is 

sum( Vi)*sum(i**4)-sum( i*i )*sum( i*i*Vi ) 


sum(l)*sum(i**4)-sum(i*i )*sum(i*i) 

where the suma are over the full range of i. For 5 points (M*2), the 
coefficients of Vi are -.0857, +.3429, +.4857, +.3429, and -.0857. For 7 
points they are -.0952, +.1429, +.2857, +.3333, +.2857, +.1429, and -.0952. 
The weighting is symmetrical about VO, the point to be corrected. The 
difference between V and VO is used as a correction factor for all range 
gates in that data scan. 
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Mot* Chat this algor itha cannot corract errors which vary with rangti 
such as thos* produced by non-horisontal scans or laser drift. These errors 
arc not asMnabla to correction. If possible they aust be recognised end 
edited out of the data prior to the application of any saoothing or oth«;r 
corrections. 



